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We discuss a state discriminator that unambiguously distinguishes between two quantum registers 
prepared with multiple copies of two unknown qubits. This device achieves the optimal performance 
by von Neumann measurement and general POVM in different ranges of the input preparation 
probabilities, respectively, and in the limit of very large program registers the optimal measurement 
reduces to only von Neumann measurements. 


The discrimination between a pair of linearly indepen¬ 
dent pure states \ipi) and 1 1 /^ 2 ) is possible only when their 
complete classical knowledge is known or these states 
carry some symmetry. Even under these conditions, their 
identification is achieved only with a certain probability 
if these states are nonorthogonal. There are strategies 
to achieve the optimal results. One is to distinguish be¬ 
tween the states with a minimal probability of error [l, 2], 
and the other one of unambiguous discrimination (UD) 
allows the possibility of a minimal inconclusive result [sj- 
[V]. Moreover, there is the minimax approach [8j to deal 
with the situations with no information about the a prior 
probabilities of the states. For practical applications, one 
needs a device that achieves the optimal discrimination 
between any pair of quantum states in a versatile way— 
first of all, it works universally on arbitrary pairs of the 
target states o be discriminated and, secondly, it should 
perform optimally in the whole range of the a prior prob¬ 
abilities or preparation probabilities possibly arising in 
the target state preparation period. 

To realize the unknown quantum states (no classical 
knowledge available) discrimination, we rely on the sym¬ 
metry from the available multiple copies of the states we 
want to identify [M3- Here we only consider the UD of 
unknown states. In Ref. [8], D’Ariano et al. propose a 
state discriminator in which the information about two 
unknown qubits is supplied in the form of a program in a 
programmble device [13|. There one copy of these qubits 
is used as the program or reference digits. This prob¬ 
lem was generalized to the situation of multiple copies 
as the reference 0G3. In Ref. [Tlj ] . by means of the 
representation theory of U(n) group, Hayashi et al. ob¬ 
tain the optimal UD solution for the averages of two uni¬ 
formly distributed unknown qudit systems prepared with 
n copies of reference plus one copy of data, while in H3 
we apply the Jordan basis method to the corresponding 
problem of two unknown qubit systems prepared with n 
copies of reference plus m copies of data which are pre¬ 
pared with arbitrary a prior probabilities. Other recent 
progress in unknown quantum state discriminations and 
their physical realization can be found in [T3 - l23j . Many 
of these studies are about the optimal UD for the av¬ 
erages of the systems. In this paper, we will construct 
a device that works on the exact input states without 
first taking the average of the inputs—the measurement 
by the device tells us which qubit it is with a non-zero 


success probability. 

We start with the generalization of the inputs in (sj- 
Here we store arbitrary n copies of qubits \ipi) and |'0 2 ) in 
the program section of the quantum registers as follows: 

I'P"’) = |V^l) 11V , 2)2 • • • |V’l)2n—l|V’2)2n|V’l)2n+l 
l^sT) = |V^l) 11^ 2)2 • • • \'lpl}2n-lH’2)2n\i>2)2n+l, (1) 

where the data, either j^i) or |?/> 2 ), at the tail is appended 
to the program as a tensor product of copies of two un¬ 
known qubits. Without loss of generality we place 
in the odd positions and |?/; 2 ) in the even positions of the 
inputs. We need to find an optimal measurement to un¬ 
ambiguously distinguish between these inputs, bearing in 
mind the data qubits, 

1 4>i) = cos(0i/2)|O) + s\n.(0i/2)e l ^ i \l) , (2) 

where i = 1,2, are randomly distributed on Bloch sphere 
with no available information about Of s and </Vs. Gener¬ 
ally these quantum registers are assumed to be prepared 
with a prior probabilities rji and ?y 2 , respectively. 

The optimal measurement performed by our device is 
generally achieved by a positive operator value measure 
(POVM), and we denote its element of unambiguosly de¬ 
tecting I’F™) as Hi, that of unambiguosly detecting |H/| n ) 
as n 2 , and that corresponding to failure as n 0 . They sat¬ 
isfy the unit decomposition: 

i = ni + n 2 + no. (3) 

The probabilities of successfully identifying these two 
possible input states are 

(M-rinrl(^r|H2|^)=p 2 , (4) 

and the condition of no error implies that 

n 2 |^i n ) = o ni|^ n ) = o. (5) 

Without any knowledge about \ifi) and |t/> 2 ) them¬ 
selves, what we can use is their symmetry, i.e. lU/™} 
is invariant under the action of symmetric group on the 
odd positions and the tail position, while IT™) invariant 
under the corresponding action on the even positions and 
the tail. To construct Hi and n 2 that perform the UD of 
the inputs, we define the following orthornormal basis: 
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where C k is the number of ways to choose k objects from 
a group of n objects without regard to order. In terms 
of these basis, the tensor product of n copies of a qubit 
is expanded as follows: 


probability of success detection as 

Pi = <*rinii*r> = ci - ciW n \p R ,T ® i s \*r) , w 


(cos(0/2)|O> + sin(0/2)e+l))® n 

= ^cos"- fc (fl/2)sin fc (fl/2)e < ^ A /c*|e fc ). (7) 
fc=o 


Thus we construct the following POVM elements satis¬ 
fying the unambiguity (no error) in measurement: 


III = ci(Ie,t — Pe.t) ® lo 

n 2 = C 2 (/o,T — Po,T ) ® Ie, (8) 


n+1 

where P e ,t = I] | e k ) E ,T E,r{e k \ (resp. P 0 ,t = 
k =0 


n+1 

| e k)o,T 0,7++ is the projection operator onto the 

k —0 

totally symmetric subspace with respect to the n even 
(resp. odd) positions and the tail position, Ie,t (resp. 
Io,t ) the direct product of unit matrices in these posi¬ 
tions, and Ci, c 2 the non-negative parameters arising from 
the unequal a priori probabilities of +i") and + 2 "). Io 
and I E are the direct products of the unit matrices in 
the odd and the even positions of the input states, re¬ 
spectively. 

Using these elements + , where i = 1,2, we find the 


where R = E,0 but S = 0,E. We apply two tricks 
to calculate the average of the projection operators more 
conveniently. First, by the expansion 


/ C k / C k ~ 1 

I efc+,T = W^+-|e fc ) fl |0) T +—|e fc _i + |l) T , 

( 10 ) 

for k = 1, 2, • • •, n, we rewrite the projection operators as 


Pr,t ® Is 

= {|eo)fl|0)T t+o| r(0| 

n I ni / rii- 1 

+ E+ T+H + k+t + 1/7 T —|e;-i)fi|l)r) 

1=1 V ° n + 1 V U n+1 

x cT - R ^ + \J 1 

+ |en)fl|l)'T R( e n\ t(1|}®+. 

(11) 

Second, we apply Eq. (7) to expand the input states into 
two parts: 


+r> = E cos"-^/2) sin++2)e^ (cos(0 i /2)|O) T + sin(^/2)e^|l) T ) 

1—0 

--- > 

i*;> 

x Ecos' l - fe (0 i /2)sin fc +/2)e ife ^ v / E|e fe ) 5 , 

k—0 


( 12 ) 
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where i = 1,2 but j = 2,1. Putting the expectation 
values from these parts together, we obtain 


(*r\p R ,T®is\*r) 

n — k + 1 


E 

k -0 


n +1 


cos 2 fa/2) 


k + 1 
n +1 


sir? (Oil 2) 


x cos 2 < n - fc) (0j/2) sin 2k (0j/2) 

” 0 7. 

+ E — I C'*c O s 2 (“- fc )(0 j /2)«n 2fc (^/2)c O t(0 i /2) 

fe=i 

x cosfa/2) sin(9i/2)cos(<j)i — 4>j) . (13) 


The total success probability P to discriminate between 
the two unknown states, assuming |\Efa) is produced with 
a priori probability ??i and fa™) with 772 , is given by 


p = mpi + r i2P2 ■ 


(14) 


Since we have no knowledge about the data states, the 
performance of the device performs is indicated by the 
average 


1 ^ r27V f' 7T 

P = J Q d( t>i J Q dOjsinOMmn^T) 

+ V2(n n m% n )) 


Tl 

= fa ci +P 2 C 2 ) . (15) 

2 (n + 1) 

of the success probability. We want to maximize this 
expression subject to the constraint that IIo is a positive 
operator. 

Let H be the Hilbert space of the input 
states given in Eq.(l), which, as it is seen 
from Eq.(lO), is spanned by the orthornormal 
basis {\ei) o \e m )E\0)T, \ei)o\e m ) e\?)t}, where 
l,m = 0, l,---,n. To have a nice matrix represen¬ 
tation of the inconclusive operator, 

n 0 = / - ni - n 2 

= (1 — Ci — c 2 )/ + c\Pe,t ® lo + c 2 Po,t ® Ie (16) 
we apply the following orthogonal transformations 
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-|ei)o|e m -i)E|l)T, 


n+1 


(17) 


where l = 0,1, • • •, n, and m = 1,2, • • •, n. The trans¬ 
formed vectors fa m ) and \xim) satisfy (?? i m \xi'm ') = 0, 
(jllm\r}l' m') — and {XlmlXl'm') = 

If we put them together with the untransformed 
|ei)o|eo)s|0)T and |e;)o|e n )£;|l)T, for l = 0, l,---,n, 
all these unit vectors will form the orthornormal ba¬ 
sis, {\rjim), \Xim), |ei)o|e 0 )_E|0) T , |ei)o|e„)_E|l)T}, of H. 
With this basis, the operator n 0 is represented by the 
direct sum of two series of block diagonal matrices: 


n 0 = n< 1} © n^ 2) , 


where 
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(i) 

0 
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n 


( 2 ) 
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V JnJ 

( K ° \ 
fa 
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V K n J 
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(19) 


These block diagonal matrices are generated regularly; 
Jo and Kq are lxl matrix 1, and J\ is a 3 x 3 matrix: 


Ji 


/ 1 1 V5T 

n+1 ^ (n+1) 3 / 2 ^2 

Vn 1 n „ 

(n+1)3/2 C 2 1 (n+1) 2 C 2 

_ n n 3/2 

\ („+l)3/2 C 2 („+l)2C 2 


_ n 

(n+1) 3 / 2 C2 

n 3 / 2 

I^+TF C2 

1 n 2 

1 ~ Cl (1+nj 2 C 2 


\ 

/ 


(20) 


I Vim) 



|ez)o|e m )s|0) 


T 


and then the size of them grows up with the general J/ 
and Ki given as (21 + 1) x (21 + 1) matrices: 
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(n+1) 3 / 2 C2 


v^TT^I) 

(n+1) 3 / 2 C2 

l + T^TFC 2 
(n+2— 2l)y/n ^ 
(n+1) 2 C2 


1 - Ci 


^//n(n+l — Z) 
(n+1) 3 / 2 C2 
(n+2— 2l)y/n 
(n+1) 2 C2 
, (i— l)n+l—Z—n 2 
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0 

0 
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1 „ I ln—n—n 2 „ I 

1 - Cl + (ra+1) 2 C 2 j 


( 21 ) 


A'l = 


1 n+l C 2 

-^//(n+l —Z) 
( ra +l)3/2~C 2 
-y/ Zn(n+1 — Z) 
(n+1) 3 / 2 C2 


\/ Z(n+1 —Z) 
(n+1) 3 / 2 C2 


1 1 l—nl—1 „ 

1 + (n+1) 2 C 2 

(n+2—2Z)+n 
(n+1) 2 


C 2 1 — Cl + 


yj Zn(n+1 —Z) 
(n+1) 3 / 2 C 2 
(n+2— 2l)y/n 
(n+1) 2 C2 
(Z— l)n+l— l—n z 
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0 

0 

c 2 • • • 0 


1 „ 1 ln—n—n 2 „ / 

' ' ' 1 - Cl + ( „ +1) 2 C 2 ) 


( 22 ) 


r 


All these real symmetric matrices can be diagonalized. 
To guarantee the non-negativity of n 0 , we only need to 
let all of their eigenvalues, which are real due to the sym¬ 


metry in the above matrices, be 11011 -negative. Through 
the induction on n, we obtained the eigenvalues of J; or 
K{ in the diagonal elements of the following: 


/ 2 -C 1 -C 2 

2 


J'l = K\ = 


\ 



(n 2 —2n—l)ciC 2 
2(n+l) 2 


0 


2 —Cl—c 2 
2 



(n 2 —2n—l)ciC 2 
2(n+l) 2 


0 


0 


0 N 
0 

1 / 


(23) 


The constant eigenvalue 1 exists for all J/’s and K{ s, 
where l = 0,1 , All the rest eigenvalues exist in 

pairs, and the sum of each pair is 2 — ci — c 2 . After ar¬ 
ranging the obtained eigenvalue pairs in all J(’s and K ['s 
by the members with the minus sign before the square 
root in an ascending order, we see that the least pair, 

1 ~ I ci “ 5 C 2 ± yj |c? + + ra 2 2 (n++ ~ c i c 2» are the com " 

mon eigenvalues of <// and Ki for l = 1 , 2 , • • •, n. and the 
second least, 1 - \c\ - |c 2 ± ^+1 + +| + n 2 (n+"]~ 2 -cic 2 , 
the common eigenvalues of J\ and Ki for l = 2,3, • • • .n, 
and so on. Therefore, if we require that the least com¬ 
mon eigenvalue of all matrices except for Jq and Kq (the 
eigenvalue of which is the constant 1 ) be non-negative, 
i.e. 


non-negative. 

Since (0111^|<^>) < 1 for V|</>) € H and all n, we have Ci < 
1 for i = 1,2 and a non-negative number 1 — +1 — |c 2 . 
Then, from the above inequality, we express c 2 in terms 
of Ci: 


C2 < 


1 - Cl 


1 - 


2n+l 

(n+1) 2 


Cl 


(25) 


To achieve the maximum success probability, we choose 
the equal sign, i.e. let the least eigenvalue of IIo be 0 . 
Inserting the resulting expression into Eq. (15) gives 


P = (?7iCi + m 


1 — Ci . n 

1 - (f+^ci 2+TI) ' 


(26) 



„ 1 , n 2 — 2 n — 1 ^ „ 

C ‘ + ^ + 2 ++ 1 P ClC2 - ° 


the inconclusive operator IIo will be guaranteed to be 


We can easily find Ci )0p t where the right-hand side of this 
expression is maximized and the corresponding c 2 j 0 p t: 


^1, opt 


(n + 1) 2 n j I — 771 

2n+l 1 n +IV Vi ' 
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FIG. 1: Optimal average success probability, P, vs. the a 
priori probability, 771 , for n = 2. Dashed line: Pi from Eq. 
(31), dotted line: P 2 from Eq. (32), and solid curve: Ppovm 
from Eq. (28). The optimal P is given by P 2 for 77 1 < 4/13, 
by Ppovm for 4/13 < 771 < 9/13 and by Pi for 9/13 < 771 . 
The lower bound if the optimal average success probability 
reaches 0 . 2 . 


FIG. 2: Optimal average success probability, P, vs. the a 
priori probability, 771 , for n = 6 . The optimal P is given by 
P 2 (dotted line) for 771 < 36/85, by Ppovm (solid curve) for 
36/85 < 771 < 49/85 and by Pi (dashed line) for 49/85 < 771 . 
The lower bound if the optimal average success probability 
reaches 0.23. 


C 2 ,opt 


(n + l) 2 
2n + 1 


( 1 - 


Vi 


1 V 1 - m 


( 27 ) 


Substituting these optimum values into (15) gives the 
optimum result for POVM: 

PpovM{n, 771) = 2 (n + 1 - 277^771(1 - 771)) . ( 28 ) 

Then we have the following optimum POVM elements: 


measurement, which projects the input onto the compli¬ 
ment of the totally symmetric subspace with respect to 
all the even digits in the program and the data digit. 
Then, we get the success probability for each operator: 
Pi,opt = m(lE,T ~ Pe,t ) ® Io\W) and P 2 ,o P t = 0, in¬ 
dicating the second input is sacrificed completely. These 
results yield the average success probability, 

_ fl 

Pl( ".’ n > = Ti 2(^TT)’ (31) 


ni, opt (n,77i) = 


x 

n 2 ,opt(n, 771) = 

x 


(n + l) 2 n / 1 — 771 

2n+l n + 1 V Vi 
( Ie,t — Pe,t ) ® lo 

( n + 1 ) 2 (l n I Vi \ 
2 n+l ( n + 1 V 1 - Vi 

( Io,t — Po,t ) O Ie ■ ( 29 ) 


Here we have assumed that a prior probabilities of the in¬ 
puts satisfy 771 +772 = 1, i.e. no failure in the preparation 
period. 

This optimal POVM, however, holds valid only within 
the range of p\ or p 2 where both 0 < Ci j0pt < 1 and 
0 < C 2 , op t < 1 are satisfied simultaneously. From Eq. 
(27) we can find the range: 


i 2 + (n + l) 2 


< 771,772 < 


(n + l) 2 


n 2 + {n + l) 2 


( 30 ) 


^ + ( ??2 = n‘ 2 + (n+ 1) 2 )> Hi,opt = ( d E,T — 

Pe,t)®Io and n 2j0 pt = 0 . The continuity to the outside 
of the POVM’s validity domain requires this structure 
remain for 1 > 771 > ■ In other words, when 

the preparation is dominated by the first input, the opti¬ 
mal measurement is realized by standard von Neumann 


for pi > ■ Conversely, for 77 ! = n , + ^ +l) , (p 2 = 

n^+^n+i) 2 ')' n 2 ,opt = ( Io,t — Po,t) ® Ie and U hopt = 0 . 

Also from the continuity, we require this structure remain 
2 

for 0 < 771 < n ‘i + ^ n+ iyi ■ This is a standard von Neumann 
measurement projecting the input onto the compliment 
of the totally symmetric subspace with respect to all the 
odd digits in the program and the data digit. Since the 
first input is sacrified completely in this case, we have 


P 2 {,n,p 1) = p 2 


2(n + l) 


= (1 - Vi) 


2 (n + 1) ’ 


(32) 


2 

for 0 < 771 < w 2 + ”„ +1 ^ • The situation is fully symmetric 
in the a priori probabilities of inputs. In the intermediate 
range, where neither of the input states dominates, the 
POVM does the job of discrimination between them bet¬ 
ter than von Neumann measurement, while outside the 
range two full projectors do the work of identifying them 
at best. The optimal average success probability of the 
device working in the whole range of a priori probability, 
as the generalization of the results in |9J, is summarized 
as follows: 

P° Pt {n,pi) 
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f P 2 (n,?7i) 

< PpovM(n,r] i) 

[ Pi (n, 771) 


® < nT+fra+l ) 7 > 

» 1 2 3 4 < _, < ( rt + 1 ) 2 
nHIn+l ) 2 — '! — n 2 + (n+l ) 2 ’ 

(ra+ 1) 2 . < 

n 2 + («+l) 2 < ^ - X - 


(33) 


It is seen from the above results that the optimal success 
probability changes continuously at the boundaries of the 
respective regions of validity for the three measurements. 
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FIG. 3: Optimal average success probability, P, vs. the a 
priori probability, 771 , for n = 00 . In this upper bound limit, 
the optimal scheme reduces to two full projector probabilities 
Pi (dashed line) and P 2 (dotted line). The minimum value 
can reach 0.25. 

Some features should be noticed for this device: First, 
from Eq. (30), the validity range of the POVM mea¬ 
surement will become narrower and narrower with the 
increase of n, the number of copies used in the program. 
As n tends to infinity, this range will shrink to a point 
at 771 = 772 = 0.5. It means that the optimum measure¬ 
ment will reduce to only von Neumann in this extreme. 
Second, if the preparation probability is fixed, the suc¬ 
cess probability of the device will increase with the num¬ 
ber of copies of the states stored in the program. The 
larger size of the program yields more information about 
the unknown states. For example, at 771 = 772 = 0.5, 
where the identification of the input state is the hard¬ 
est and the difference between the POVM and two full 
projectors is the largest, the average success probability, 
PpovM(n) = n/(4n + 2), increases from 1/6 to 1/4 by 
50% as n goes from 1 to infinity. We depict these features 



in Figs. 1-3 for three different n {n = 2, 6, and 00 ). 

Another interesting feature is that the optimal mea¬ 
surement operators and the validity domain depend only 
on the number of the states we put into the program 
and their preparation probabilities. Since they are inde¬ 
pendent of the information encoded in a specific pair of 
unknown qubits |t/>i) and ^ 2 ), the device performs uni¬ 
versally. It can be programed with whatever couple of 
unknown qubits we want to identify and discriminates 
between the produced registers with the best chance of 
success and no error. 

Finally, we compare the obtained results with those 
of UD for the averages of systems. Suppose the a prior 
probabilities of the systems prepared as in Eq. (1) are 
equal. If we take the averages of the density matrices 
Pl = and p 2 = I'&rX*!"! over the Bloch 

spheres of the uniformly and independently distributed 
1 7 / 1 ), 1 7 / 2 ) and find the optimal UD solution to these 
two averages, which are two uniformly distributed known 
mixed states, we will have the optimal average success 
probability of 1/3 if the number of the qubit copies tends 
to infinity (TT], (TJ . This is equal to the average of IDP 
limit [3]-[5] for discriminating a pair of known qubits. 
The corresponding maximum average success probabil¬ 
ity, 1/4, of our device is obtained through first measur¬ 
ing the systems with the optimal POVM and then tak¬ 
ing the averages of the expectation values of III and II 2 , 
and is lower than this IDP average. However, the device 
described in the current work realizes a practical iden¬ 
tification of any pair of randomly distributed unknown 
qubits instead of their averaged input states. We see 
from this fact that the discrimination of two unknown 
systems (with only some information about their sym¬ 
metry available) and the discrimination of their known 
averages are not equivalent. 

So far we have worked out a scheme for a device that 
unambiguously discriminate between pairs of quantum 
registers produced with more than one copy of unknown 
qubits used in the program or reference digits. Its maxi¬ 
mum average success probability increases with the num¬ 
ber of the copies we store in the program, and also de¬ 
pends on the a prior probability of the states to be dis¬ 
criminated. We can tune the parameters ci and C 2 of the 
setup according to the copy number n and the a prior 
probability rji or 772 , so that it optimally performs in the 
required measurement. 

The authors would like to thank M. Hillery for valuable 
discussions. 
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